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ABSTRACT 

Statistical observations of the epoch of reionization (EOR) power spectrum provide a rich data set for under- 
standing the transition from the cosmic "dark ages" to the ionized universe we see today. EOR observations 
have become an active area of experimental cosmology, and three first generation observatories — MWA, PAST, 
and LOFAR — are currently under development. In this paper we provide the first quantitative calculation of 
the three dimensional power spectrum sensitivity, incorporating the design parameters of a planned array. This 
calculation is then used to explore the constraints these first generation observations can place on the EOR 
power spectrum. The results demonstrate the potential of upcoming power spectrum observations to constrain 
theories of structure formation and reionization. 

Subject headings: Cosmology: Early Universe, Galaxies: Intergalactic Medium, Radio Lines: General, Tech- 
niques: Interferometric 
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1. INTRODUCTION 

The reionization history of the universe provides an im- 
portant tool for understanding the epoch of structure forma- 
tion and the appearance of the first luminous objects. Ex- 
isting experimental clues about this period are confusing, 
however. The large Thompson scattering measured by the 
WMAP satel lite implies that reionization occurred by red- 
shift z ~ 15 (Kogutetal. 2003), whereas quasar absorption 
line st udies show significant neutral hydrogen a t redshift 
6 JDiorgovski et alJ 120011 iBecker et alJlfooH iFan et all 



2003; Wvithe & Loeb 2004b). Reconciling these measure- 
ments s eems to require a fairly complicated reionization 
history (lHaiman & HoldeJ 120031 ICenl 120031 ISokasian et aTJ 
120031 iMadau et all2004l) 

Directly observing the process of structure formation dur- 
ing the epoch of reionization (EOR) would resolve these un- 
certainties. As the primordial hydrogen cools and later re- 
heats, density contrasts are revealed as fluctuations in the 
brightness tem perature of the re dshifted 21 cm neutral hydro- 
gen line ijS unvaev & Zeldovi chlll972l iH ogan & Reeslll979t 
I Scott & Reeslll990t [S hapiro et al. 1994: lLoeb & Zaldarriaga 
120041 iBarkana & Loebl l2005bl) . Additionally, as the first 
luminous objects ionize their surrou ndings bubbles ap- 
pear in th e diffuse 2 1 cm emission ( Madau et alJ 119971 
[ Tozzi et alJ 120001 ICiardi & Madaul 120031 IZaldarria ga et all 
l2004UFurlanetto et all2004al) . 

The neutral hydrogen emission from this period appears 
in the low radio frequencies as faint fluctuations between 
75 — > 200 MHz (for redshifts 18 — > 6). Directly imaging the 
fluctuations will require the sensitivity of the Square Kilo- 
meter Array ( F urlanetto & Briggsl2004l) . but statistical obser- 
vations of the fluctuation power spectrum should be obtain- 
able with the much smaller first generation EOR observato- 
ries currently under construction. Measuring the power spec- 
trum and its evolution would provide a wealth of information 
about structure formation and the fundamental astrophysics 
behind reionization. At high redshift, observations would ex- 



plore structure formation in the line ar gravity regime and con- 
strain the properties of dark matter ( Loeb & Zaldarriaga 2004; 
iPenl 12004). and at lower redshifts they would probe reion- 
ization and follow the emerge nce an d propertie s of the first 
luminous objects jTozzi e t al.ll200 0l ISomerville et al.l 120031 
iFurlanetto et all2004albUlliev et all2005l) . 

Measuring the EOR power spectrum builds on the sta- 
tistical techniques developed for analyzing the cosmic mi- 
crowave background (CMB) anisotropy. Unlike the CMB, 
however, the EOR signal is fully three-dimensional since the 
frequency of the redshifted 21 cm line maps to the line-of- 
sight distance. Recent efforts have shown how the three- 
dimensional nature of the EOR signal can be used to increase 
instrumenta l sensitivity an d mitigate against foreground con- 
tamination Morales & Hewitdl2004l IZaldarriaga et al.ll2004l 
iBharadwai & Alil2005l) .~ 

In th is p aper, we utilize t he formalism of Moral es~& Hewitt! 
J2004I) and lMoralesI J2005I) with the parameters of a planned 
array to provide the first quantitative calculation of the three 
dimensional power spectrum sensitivity of forthcoming ob- 
servations. The calculations include realistic observational 
parameters such as antenna layout, field of view, and antenna 
temperature, and provide a fiducial mark for the capabilities 
of the first generation EOR observations. 

We begin in Section 2 by detailing the observational param- 
eters and techniques used for our calculation. Section 3 then 
presents the results, analyzing the EOR sensitivity as a func- 
tion of length scale, redshift, and global ionization fraction. 
We conclude in Section 4 with a discussion of the additional 
experimental issues which must be confronted in a realistic 
array. 

2. FIRST GENERATION EOR EXPERIMENTS 



The first generation of radio observatories targeting the 
EOR power spectrum consists of three instruments that are 
currently under development and should be operational by the 
end of the decade. The instruments are PAST, MWA, and LO- 
FAR 

Electronic address: jdbowman@mit.edu; mmorales@space.mit.edu; jhewitt@space.mitjqdjj p r i meV al Structure Telescope (PAST) is a project of 
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FIG. 1. — This photograph shows the first MWA prototype antenna being 
tested on site in Mileura. The antenna consists of 16 crossed-dipoles in a four- 
by-four meter grid and is elevated approximately 0.5 m above the ground. 
Five hundred antennas will be placed in a 1500 m diameter array for the low 
frequency band of the first-generation MWA. 



TABLE 1 

Fiducial Observation Parameters 



Parameter 


Value 


Array configuration, p(r) 


~r- 2 


Array diameter, D 


1500 m 


Bandwidth, B 


8 MHz 


Frequency resolution 


8 kHz 


Integration time, t 


360 hours 


Number of antennas, N 


500 



NOTE. — Observation parameters used in the sensitivity analysis. 



TABLE 2 
Redshift Dependent Parameters 



polarization yagis observing the north celestial pole. PAST 
will b e the first EOR array to begin acquiring data. tPen et all 
2005, submitted) 

The Mileura Wide-field Array (MWA) is a collabora- 
tion of the Massachusetts Institute of Technology, Harvard- 
Smithsonian Center for Astrophysics, the Australian Univer- 
sity Consortium, the Australian National Telescope Facility, 
and the Western Australian government. The MWA will be lo- 
cated at the Mileura radio quiet site in the Western Australian 
desert. The low frequency portion of the array will consist of 
8,000 dual-polarization dipoles and feature very high spatial 
dynam ic range and c alibration precision, including polariza- 
tion. (Salah et al. 2005, in press) 

The Low Frequency Array (LOFAR) is being built by AS- 
TRON in the Netherlands, and will consist of 72,000 dual po- 
larization dipoles operating at EOR frequencies. While the 
radio interference environment of the Netherlands is challeng- 
ing, LOFAR will have by far t he largest collectin g area of the 
first generation EOR arrays. (Kass im et al.l20 04) 

2. 1 . The Fiducial Observation 

To calculate accurately the sensitivity of an EOR measure- 
ment, we need to specify both the details of the instrument 
and the observing strategy. For this paper we wish to choose 
a set of realistic observation parameters so that we can discuss 
the statistical sensitivity of first-generation EOR experiments. 

We have chosen our reference measurement to be a deep 
observation of a single target field using an array configura- 
tion based on the MWA. Of the first generation arrays, the 
MWA is intermediate in theoretical sensitivity and has good 
control of systematic errors. 

Many factors contribute to the instrumental response of an 
array. Among the most important to consider are the field of 
view, angular resolution, collecting area, antenna distribution, 
and bandwidth. 

For the MWA, each antenna consists of 16 crossed-dipoles 
in a four-by-four meter grid as shown in Figure ^ A beam- 
forming operation on the dipoles within an antenna produces 
a field of view of 15-45° (FWHM), depending on frequency. 

The reference array consists ofN = 500 antennas distributed 
within a D = 1500 m diameter circle. The density of antennas 
as a function of radius is taken to be approximately p(r) ~ r~ 2 . 
The angular resolution is given by X/D and the collecting area 
by N dA, where dA is the collecting area of each antenna and 
scales like dA ~ 16(A 2 /4) for wavelengths below two meters. 

Although the bandwidth of the MWA is 32 MHz, we re- 





z = 6 


Z = 8 


z= 10 


Z=12 


Angular resolution (°) 


0.06 


0.07 


0.09 


0.10 


Antenna collecting area, dA (m 2 ) 


9 


14 


18 


18 


Antenna response scale, (°) 


19 


31 


38 


43 


Frequency (MHz) 


203 


158 


129 


109 


System temperature, T„ s (K) 


250 


440 


690 


1000 



NOTE. — Characteristics of the fiducial observation that depend 
on frequency, and thus on redshift. Note that the antenna collecting 
area is capped above z = 10 due to self-shadowing by the antennas at 
low frequencies. Also, the system temperature is dominated by sky 
noise and depends significantly on frequency. 



strict the bandwidth of our reference observation to B = 8 
MHz. This avoids complications introduced by cosmic evolu- 
tion, yet still provides measurements along the line-of-sight 
at the length scales of the strongest fluctuations in the 21 
cm signal. It is also th e minimum bandwidth suggested by 
Wvithe & Loeb (2004a) to ensure sensitivity to fluctuations 
during the final stage of reionization when ionized bubbles 
first overlap completely. 

Finally, the integration time on the target field must also be 
specified. Observing for four hours per day over the course of 
a six month season would result in approximately 720 hours 
of integration. Conservatively rejecting half the data for non- 
ideal conditions yields 360 hours of integration during the 
most favorable circumstances. 

The full parameter set for our fiducial observation is sum- 
marized in Tables [Qand|2] 

2.2. The Data Cylinder 

Since neutral hydrogen is optically thin to the 21 cm line, 
the visibility measurements of EOR observatories inherently 
sample the emission from a three dimensional volume of 
space at high redshift. To good approximation, these mea- 
surements form a three dimensional data cylinder in visibility 
space (u,v,f) due to the overall circular shape of the array. 

By applying Fourier transforms along one or more of the 
coordinates of the data cylinder, the measurements also may 
be represented as cylinders in several additional useful coor- 
dinate spaces. These include real space (with units of comov- 
ing Mpc), cosmological Fourier space (k = k\,k2,ki), image 
space (8 x ,9 y ,f), and instrumental Fourier space (u, v, if). Each 
coordinate space possesses advantages for different stages of 
the analysis. For example, point source foreground removal 
is most conveniently accomplished in the image space, while 
the power spectrum of the EOR signal has symmetries that 
are most easily exploited in the Fourier space. Since space 
is isotropic (rotationally invariant) the EOR signal is approx- 
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TABLE 3 

Data Cylinder Dimensions for Redshift z = 8 



Frame 


Diameter 


Depth (Line-of-sight) 


Cosmological Fourier (k) 


0.55 Mpc" 1 


48 Mpc -1 


Image (0 X , G y , /) 


62° 


8 MHz 


Instrumental Fourier (n,v, rf) 


790 A 


1.2- 1CT 4 Hz" 1 


Real 


2400 Mpc 


130 Mpc 


Visibility (u,v,f) 


790 A 


8 MHz 



NOTE. — Dimensions of the data cylinder at redshift z = 8 for 
the example array configuration in different frames. The diameter of 
the cylinder is in the plane of the first two coordinates and the depth 
is along the third coordinate. 



imately spherically symmetric in the Fourier space. All of 
the measurements from a spherical shell are thus drawn from 
the same statistical ensemble and can be averaged to gether 
to maximize t he signal to noise tMorales & HewitJ 120041: 
Morales 2005). The spherical symmetry is the basis of the 
three dimensional statistical EOR measurement and employed 
in our sensitivity calculation below. 

At redshift z = 8, the data cylinder in cosmological Fourier 
space has a diameter of 0.55 Mpc" 1 in the fcjA^-plane and 
spans 48 Mpc -1 in the line-of-sight direction. Clearly the data 
cylinder is very elongated in Fourier space, with the line-of- 
sight axis responsible for the high spatial frequency contri- 
bution. Table [3] lists the dimensions of the data cylinder in 
additional frames. We will see in the following sections that 
the dimensions of the data cylinder play an important role in 
the planned statistical measurements. 

2.3. The Instrumental Response 

We have so far outlined the details of the fiducial observa- 
tion. In order to understand the sensitivity of the experiment, 
the 21 cm power spectrum must be mappe d to an instr umen- 
tal resp onse. Following the development of Morales & Hewitt 
(2004|, this response is given by the convolution of the power 
spectrum, P(k), with the instrumental window function, W(k), 
according to 

C 7 (k) = (|A/(k)| 2 )= J P(k)\W(k-k')\ 2 d 3 k'. (1) 

The window function is given by the array's spatial and fre- 
quency response. We approximate it with a function that is 
independent of frequency within the band and that depends 
on 9 as 

w(<9) = cos 2 (§6ye), e <e (2) 

where 9 is the antenna response scale size and W(9) is zero 
outside the defined region. In general the Fourier space win- 
dow function, W(k), is a very sharply peaked function. At 
redshift z = 8, the width of W{k) for the example array config- 
uration is approximately 0.003 Mpc" 1 in the fc^-plane and 
0.05 Mpc -1 along the ^-direction. 

Inherent in every radio observation is thermal noise. The 
contribution of this noise can be estimated by dividing the 
Fourier space into a large number of independent cells, where 
the size of each cell is given by the window function, W(k). 
The thermal noise per indep endent cell can then be approxi- 
mated using (Morales 2005, His Equation 11) 
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FIG. 2. — Four model antenna configurations (top), the corresponding den- 
sities of visibility measurements (middle), and 1-cr uncertainties in the mea- 
sured power spectrum due to thermal noise (bottom). The antenna configura- 
tions are characterized by power-law density profiles, p(r) ~ r~ 3 (solid), r~ 2 
(dash), r~' (dot), and r~° (dash-dot). The abscissas of the three panels are 
aligned so that the r, u, and k coordinates correspond. The vertical gray bars 
represent the 4 m width of an antenna (far left), the 750 m radius of the ar- 
ray (middle) and the 1500 m maximum baseline in the nv-plane (right). The 
uncertainties in the bottom panel extend beyond the bounds of the maximum 
baseline due to the elongated k} axis, as described in Section l2~2l 



where dA is the physical antenna area, drj is the inverse of the 
total bandwidth, Icb is Boltzmann's constant, T sys is the total 
system temperature, e is the efficiency, B is the total band- 
width, n is the time average number of baselines in an observ- 
ing cell, u = V u 2 + v 2 , and t is the total observation time. 

The thermal noise is approximately independent of fre- 
quency within the observing band and has a cylindrical sym- 
metry in Fourier space. The measured value of the power 
spectrum, however, is an average over spherical shells in 
Fourier space. Thus, to calculate the uncertainty due to ther- 
mal noise in the measurement, we must average the uncer- 
tainty contributions from all of the independent cells in spher- 
ical shells. The result is that the uncertainty at a given length 
scale has a somewhat complicated dependence on thermal 
noise, and therefore, on the antenna distribution of the array. 

The top panel of Figure [2] shows four candidate antenna 
profiles for the MWA. All of the distributions contain 500 
antennas and respect the minimum physically possible an- 
tenna spacing (which is responsible for the density plateaus 
as r — ► 0). The four distributions are labelled by the antenna 
density profiles at large radii: p(r) ~ r°, r , r" 2 , and r" 3 . 

The second panel shows the density of visibility measure- 
ments in the Mv-plane for each of the antenna distributions. 
The density is related to the average number of measurements 
per cell, n(u), in Equation[3]bythe cell size. 

The third panel of Figure |2] gives the uncertainty due to 
thermal noise per spherical logarithmic shell in cosmological 
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Fourier space with five bins per decade. Two limiting regimes 
are evident in this plot. For small shells in Fourier space, the 
noise depends strongly on the density profile since only visi- 
bility measurements from baselines smaller than a shell's ra- 
dius contribute. Increasing the steepness of the antenna pro- 
file condenses the visibility measurements toward the origin 
of the k\k2-plane, reducing the noise in small shells. In the 
other limiting case, shells in Fourier space that have extended 
beyond the radius of the largest baseline in the fc^-plane 
(k > 0.5 Mpc -1 ) include information from every antenna, and 
the density profile has little effect on the noise. 

Second order effects such as the distinction between coher- 
ent integration of visibility measurements within a cell and the 
incoherent averaging of independent cells leads to the small 
differences seen at large k. Between these limiting cases, the 
interaction of the cylindrical symmetry of the thermal noise 
with the spherical symmetry and logarithmic widths of the 
shells causes a more complicated behavior. At length scales 
comparable to the diameter of the data cylinder, where visibil- 
ity measurements are sparse over much of a shell, the uncer- 
tainty even has a local maximum for the centrally condensed 
arrays. 

It is clear from Figure |2] that the sensitivity of the statisti- 
cal EOR measurement is tightly linked to the array antenna 
distribution. The difference in uncertainty between the uni- 
form distribution and the steeper power-law distributions is as 
much as two orders of magnitude, depending on the length 
scale of interest. As indicated in Table ^ we will use the 
p(r) ~ r~ 2 antenna distribution as our reference in the remain- 
der of the paper. 

2.4. The Model Power Spectrum 

For our calculations we use a simple model of the redshifted 
21 cm power spectrum that has been used commonly in the 
literature, in which the hydrogen in the intergalactic medium 
(IGM) is fully neutral, follows the dark matter distribution, 
and ha s a spin temperature much larger than the CMB temper- 
ature (IMaSiZiEiDlillll^iilillulMoO; Zaldarriagaet alJ 
2004). This is a reasonable model for the fluctuations after the 
spin temperature has been heated by the Wouthuysen-Field 
effect or other processes and before reionization by the first 
luminous objects, though over what redshift range this may 
be observed is uncertain. 

This model power spectrum is computed using CMBFAST 
(Seliak & Zaldarriag alll996|) and does not include velocity 
distortions, but Barkana & Loeb (2005a) have shown that in- 
cluding peculiar velocities increases the signal amplitude by 
about a factor of two. Distortions due to geometrical ef- 
fects (such as a scaling between the line-of-sight and sky- 
plane axes) and peculiar velocities could allow separation of 
cosmological and astrophysical effects and provide sensitive 
probes of the unde rlying cosmology, but are not included in 
our simple model jAlcock & P aczvnskil ll979T: lKaiserjfi987t 
Barkana & Loeb 20QH). 

Figure [5] shows the results of performing the above calcu- 
lations for the fiducial experiment using our simple neutral 
hydrogen power spectrum at redshift z = 8 with standard con- 
cordance cosmology (fi^ = 0.3, = 0.7, h = 0.7). The power 
spectrum was convolved with the instrumental window func- 
tion to produce the solid blac k line, which is p lotted in the 
observer units Jy 2 Hz 2 (see [Morales & Hewittl (|2004) for a 
discussion of units). 

Just as there is an inherent uncertainty due to thermal noise, 
there is also an inherent cosmic sample variance in the ob- 




10 3 icr 2 10 1 io° 10 1 

k [Mpc -1 ] 

FIG. 3. — Combined 1-cr uncertainty for logarithmic shells of width k a2 
(equivalent to five spectral points per decade) in the measured power spec- 
trum due to sample variance (dark gray) and combined sample and thermal 
variance (light gray). The instrumental response to the ideal neutral hydro- 
gen power spectrum at redshift z = 8 is shown in black. Logarithmic shells 
of width k ' 2 were used for bins. The details of the example observation con- 
figuration are specified in Tables [Tl and |2| In the upper panel the signal is 
shown in observer units (where the measurement errors are Gaussian), while 
the lower panel uses the theoretical convention by changing the ordinate to 
(k^P/lir 2 ) 1 / 2 and plotting on log-log axes. 



served power spectrum. For a spherical shell in Fourier space, 
the sample variance can be estimated by considering the num- 
ber of independent samples in the shell and approximating by 
using Gaussian statistics. This uncertainty is plotted as the 
dark gray region around the ideal power spectrum in Figure 
13 while the combined uncertainty due to sample variance and 
thermal noise is added in quadrature, and the light gray region 
shows the full 1-sigma uncertainty. 



3. RESULTS 



3.1. 



The Measured Power Spectrum 

Several effects of the instrumental response are contained 
in the measured power spectrum shown in Figure[3] Since the 
shape of the power spectrum is smeared by convolution with 
the instrumental window function, both the relative amplitude 
of the peak and the distinction of the baryon bump at k ~ 
0.04 are slightly reduced. In addition, the uncertainty due to 
thermal noise and cosmic sample variance increases rapidly as 
k — » 0, creating unfavorable sensitivity on large length scales. 

The field of view is a significant factor determining the un- 
certainty at low k since the number of measurements at these 
scales is proportional to 9 2 . Since the line-of-sight depth of 
the data cylinder in real space is less than these length scales, 
the advantages of a three-dimensional data set do not apply 
at low k. Only spherical shells with radius k > 0.02 contain 
contributions from measurements with kj, 0. 

The antenna distribution of the array, as we saw in Sec- 
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tion 2, also affects the uncertainty of the measurement. Al- 
though more condensed arrays produced less uncertainty at 
large length scales (see Figure|2}, other considerations which 
rely on synthesis imaging, such as removing astrophysical 
foregrounds, are adversely affected by condensing the an- 
tenna distribution and could potentially negate any benefits 
from such a change. 

The sensitivity decreases again for large k since the power 
spectrum falls off more rapidly than the thermal noise. For the 
example array, the best range for constraining the measured 
power spectrum at redshift z = 8 is approximately 10~ 2 < k < 
10" 1 Mpc" 1 . 

3.2. Redshift Range 

We can estimate the sensitivity of the array at additional 
redshifts by modifying the parameters of the fiducial observa- 
tion. Four primary characteristics of the array are frequency 
dependent: the field of view (and thus the instrumental win- 
dow function, W), the collecting area, the angular resolution, 
and the system temperature, T sys . Table |2] lists the values of 
these parameters at four frequencies corresponding to red- 
shifts of z = 6, 8, 10, and 12. 

The changes in these parameters require both the measure- 
ment uncertainty and the instrumental response to the red- 
shifted 21 cm emission to be calculated for each redshift since 
the independent cell size, data cylinder dimensions, and char- 
acteristic thermal noise are directly affected. Figure |4] dis- 
plays the results of such calculations for observations at the 
redshifts listed in Table |2] Again, the reference signal was 
a fully neutral IGM and the measurement was averaged over 
spherical logarithmic bins of width k 2 . 

The top panel, for redshift z = 6, illustrates the measurement 
with the greatest sensitivity. Two factors contribute to this 
performance: the amplitude of the power spectrum increases 
as z — > 0, and the system temperature of the instrument, pri- 
marily responsible for the uncertainty, decreases. On the other 
hand, the field of view and collecting area are reduced consid- 
erably, thus limiting the improvement in sensitivity for lower 



redshifts. The net result is that the amplitude of the power 
spectrum increases by a factor of ~ 8 between redshifts z = 12 
and 6, while the system temperature, dominated by galactic 
foreground emission, decreases by a factor of ~ 3.5. The sen- 
sitivity is sufficiently great at redshift z = 6 that the dominant 
source of uncertainty at large scales is sample variance. Al- 
though it is unlikely that no ionization will have occurred by 
redshift z = 6, even significantly weaker signals should be de- 
tectable. 

There is a clear degradation of the measurement sensitivity 
in Figure as the redshift increases until, by redshift z = 12, 
the observation is infeasible without longer integrations or ad- 
ditional collecting area. 

3.3. Sensitivity to Reionization Models 

Based on the results of the above sections, it is not unrea- 
sonable to expect that the first generation of EOR experiments 
may be able to go beyond simple detections of high redshift 
neutral hydrogen and distinguish between different theoretical 
reionization scenarios. 

Several theoretical reionizati on models have bee n discussed 
in the literature ij Santos et alJl200 3HZa ldarriaga et all 120041: 
Furlanetto et al. 2004a b; Sant os et all2005l) . Figure |5] shows 
the results of a simulation based on the example array con- 
figura tion for r e dshift z = 8 and includes for comparison 
Furlanetto et al. (2004a b) models for several power spectra 
with different reionization fractions. The solid black line and 
shaded gray regions are the same as in Figure |4] Column 
2. It is immediately evident from the figure that the large 
changes in peak amplitude of the power spectra for these mod- 
els would be ideal for constraining the ionization fraction. 

4. EXPERIMENTAL CONSIDERATIONS 

Radio observations in the meter bands are notoriously dif- 
ficult. While this paper has focused on the theoretical sensi- 
tivity of first generation EOR observations, it is important to 
remember that there are many other observational difficulties 
which must be overcome by the experiments. In this section 
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foregrounds — which include synchrotron and free-free emis- 
sion from the Galaxy, free-free emission from elections in the 
IGM, and ext ended and point sources — were an i nsurmount- 
able obstacle SDi Matteo et"afll2002t lOh & Mackll2003l) . but 
additional studies indicate that multi-frequency observations 
and statistical techniques should provide method s to separate 
the foregrounds from the EOR sign al ( Di Matteo et alJl2004t 
Zaldarriag a et all2004: iMorales & Hewittll2004tlSantos et dl 
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Another concern is the instrumental effects related to galac- 
tic emission and imperfect instrumental calibration. Galactic 
synchrotron radiation dominates the sky at low radio frequen- 
cies and is responsible for the large system temperatures used 
in the fiducial observation. Additionally, the radiation is Fara- 
day rotated by the interstellar medium and presents a bright, 
highly structured polarized emission pattern across the sky. 
This places very tight constraints on the precision of the in- 
strumental polarization calibration. 

All of these additional considerations should be manage- 
able, but they complicate the analysis of the EOR power spec- 
trum in actual observations. After removal or correction, each 
will produce a characteristic residual statistical signature in 
the observed data cylinder. Efforts are underway to assess the 
best techniques for addressing these contaminants and to esti- 
mate their residual statistical signatures. 



FIG . 5 . — Same as Figure^ Column 2 for redshift z = 8, but the data points 
show a simulated realization of the measured power spectrum. The error bars 
are the 1— a uncertainties calculated from the thermal noise. The dashed lines 
are different values of the ionization fraction in the Furlanetto et al. 1 2004a b) 
models and are, from top to bottom, = 0.51, 0.0 (solid), 0.43, 0.38, 0.25, 
and 0.13. Note that the amplitude of the power spectrum for the fully neutral 
IGM is within those of the other ionization fractions. In general, the ampli- 
tudes of the model power spectra drop rapidly between x\ = 0.0 and 0.13, and 
then slowly increase with ionization fraction as large bubbles increase the 
contrast. 



we briefly describe some of the additional observational con- 
siderations. 

The redshifted 21 cm radiation targeted by EOR experi- 
ments falls in the frequency range commonly used for tele- 
vision, FM radio, and satellite transmission. Both PAST and 
MWA have chosen very remote locations to get away from 
the radio communications which are now ubiquitous in many 
parts of the world, and all three experiments are developing 
advanced radio frequency interference (RFI) mitigation tech- 
niques. 

Additionally, turbulence in the Earth's ionosphere refracts 
low frequency radio waves. This acts much like atmospheric 
distortions at optical wavelenghts, and must be corrected us- 
ing techniques similar to wide-field adaptive optics. 

After RFI and ionospheric distortions have been removed 
from the observations, there are still astrophysical foreground 
contaminants which are five orders of magnitude brighter 
than the EOR emission. Initial analysis suggested that these 



5. CONCLUSION 

Radio observations of the redshifted 21 cm emission during 
the epoch of reionization are capable of providing statistical 
measurements of the properties of neutral hydrogen at high 
redshift. These measurements can be carried out with the rel- 
atively small first-generation EOR instruments. Future larger 
arrays, such as the Square Kilometer Array, will be capable in 
principle of mapping individual features on the sky, and larger 
collecting areas should also allow high-precision measure- 
ments of the EOR power spectrum at higher redshifts, with 
important implications for fundamental cosmological mea- 
surements (Loeb & Zaldarriaga 2004). 

In this paper, we have illustrated explicitly the fundamental 
uncertainties for a realistic experiment. This analysis provides 
a firm foundation to motivate continued development in this 
field. 

Measurements of the reionization history are fundamental 
to understanding the evolution of the universe, and with prop- 
erly designed experiments, these measurements are feasible 
and may be obtained in the very near future. 
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work was provided by NSF grant AST-0121 164 and the MIT 
School of Science. 
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